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Abstract 

Neutrino and anti-neutrino induced single-charm production are par- 
ticularly interesting to study the strange-quark parton distribution func- 
tion and the threshold effect in the cross-section, associated with the heavy 
quark production. Over the past 30 years, many experiments have car- 
ried out these studies with complementary techniques: calorimetry, bubble 
chambers and nuclear emulsions. 

In this paper, we review these data and combine them statistically to 
extract a world averaged single-charm production cross-section for both 
neutrino and anti-neutrino. 



1 Introduction 



Neutrino and anti-neutrino induced charm-production is interesting because it 
can be used to isolate the strange-quark parton distribution function and to 
study the transition to heavy quarks. In particular, the understanding of the 
threshold behavior associated with the charm-production is critical for the ex- 
traction of sin^ 9w from neutrino deep-inelastic data. Furthermore, a better 
understanding of the inclusive charm-production cross-section is crucial for the 
background determination of future experiments aiming at the study of neutrino 
oscillations. 

This paper is organized as follow: in Section ^ we briefly discuss how charm- 
production can be studied by looking at the so-called di-lepton events and review 
all available di-lepton data. In Section ^ the inclusive charm search by using 
nuclear emulsions is discussed together with the available data. In Section Q 
the di-lepton and emulsion data are combined to compute a world averaged 
inclusive charm-production cross-section. 



2 Charm-production studies with di-leptons 

2.1 Experimental issues 

Many experiments have studied neutrino and anti-neutrino charm-production 
by looking at the presence of two oppositely charged leptons in the final state. 
This technique was firstly used in 1974 |^ when the neutrino induced charm- 
production was discovered. 

In the case of neutrino scattering, the underlying process is a neutrino 
charged-current (CC) interaction with an s or o? quark, producing a charm- 
quark that fragments into a charmed hadron. The charmed hadron may decay 
semi-leptonically producing opposite sign di-leptons through the process: 



1^1, + N — > n- + c + X (1) 

^ s + 1+ + i^i 

Analogously an anti-neutrino can interact with a s or d anti-quark, producing 
a charm anti-quark that fragments into a charmed hadron, again leading to a 
final state with two oppositely charged leptons. 

Since 1974, several experiments (||^ to have used this technique to 

study charm-production. They can be separated in two classes: calorimetric 
and bubble chamber experiments, exploiting the muonic and electronic decay 
of the charmed hadron, respectively. 

Calorimetric experiments are characterized by a massive iron target and a 
muon spectrometer to identify the muon and measure its charge. For these 
experiments pion and kaon decays constitute the main background. The high 
density of the target calorimeter minimizes this background due to the short 
interaction length of the detector. A further background reduction is obtained 
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by requiring a minimum momentum, typically > 5 GeV, for the less energetic 
muon. The drawback of such a selection is that these experiments are not able 
to search for charm-production at relatively low neutrino energies. For a typical 
calorimetric experiment it is not possible to investigate energy regions below 
15 GeV, where the slow-rescaling threshold effect is more important. 

The main characteristic of a bubble chamber filled with a mixture of heavy 
liquids (Ne-H2, freon-propane) is its high efficiency in identifying electrons. 
Therefore, they searched for charm-production by looking at fj,~e'^ events. In 
these experiments, the low threshold, Pe+ > 0.3 GeV, combined with high statis- 
tics for Ely < 30 GeV, gives good sensitivity to the slow-rescaling threshold be- 
haviour. The main background sources for these searches are 7r° Dalitz decays 
and PeCC interactions. 

2.2 Di-lepton data sample and selection criteria 

Table |l| summarises the available (anti-)neutrino di-lepton data samples. 

The data selection is similar for all calorimetric experiments, although the 
kinematical cuts can be slightly different. In the following we describe qualita- 
tively the data selection flow: 

• the event must occur in coincidence with the beam and fire the penetration 
trigger (CC interaction trigger); 

• to ensure both the longitudinal and transverse containment of the hadronic 
shower a fiducial cut is applied; 

• both muons have to be well reconstructed, with the closest approach be- 
tween them being less than a certain distance (typically 0(10 cm)). This 
cut allows to reject overlays of two CC events as well as obvious muons 
from the decay of shower hadrons. 

The di-lepton events are categorised as originating from an incident neutrino 
or anti-neutrino by assuming that the primary muon (the muon produced at the 
leptonic vertex) is the one with the largest transverse momentum with respect 
to the beam direction. 

For each event, the reconstructed muon parameters at the vertex (~??j, Ei, 
l^i — "p*j/ I Pi I, i = 1, 2) and the shower energy E^ach a-re used to compute the 
following kinematical variables: 

• ^vis = £^1 + £^2 + Ehad-, the visible neutrino energy; 

• Q\is — 2£'uis(£'i — 1^1 ■ ~i*) — rrij^, the visible negative four-momentum 
trasfer squared, where is the leading muon 3-momentum and ~T rep- 
resents a unit vector parallel to the beam direction; 

• i^vis = E2 + Efiad, the visible energy transferred to the hadronic system; 

• Xvis = Qlig/2Mvyis, the visible Bjorken x] 
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Experiment 


Technique 


I sample 


/+Z sample 


CDHS § 


Calorimetry 


9922 


2123 


LBL Coll. i 


Bubble Chamber 


0(50) 




CBNL Coll.J§ 


Bubble Chamber 


O(IOO) 




NOMAD j|] 


Calorimetry 


2714 ± 227 


ii5ii? 


CCFR [| 


Calorimetry 


4247 ± 90 


944 ± 34 


Gargamelle I M 


Bubble Chamber 


62 ± 10 




Gargamelle II j8| 


Bubble Chamber 


0(50) 




E53A+E53B Ji 


Bubble Chamber 


602 ± 208 




Foudas et al. Jld] 


Bubble Chamber 


1460.4 ±42.1 


223.5 ± 5.0 


E616+E701 Inf 


Calorimetry 


852 ± 77 


68± 16 



Table 1: Available (anti-)neutrino di-lepton data samples. The samples are 
background subtracted. 

/-E^jg, the visible Bjorken y. 

In addition to the topological cuts listed above, a set of kinematical cuts 
is also applied. Both muons are required to have an energy greater than 4- 
10 GeV, primarly to reduce the meson decay background. The visible hadronic 
energy is required to be more than 5-20 GeV to ensure a good hadronic energy 
reconstruction. 

The bubble chamber experimental search for charm-production relies upon 
its high efficiency in identifying electrons. Therefore, in order to identify /i^e^ 
dilepton candidates, the film is scanned for all events with an e"*" coming from 
the vertex with a momentum greater than 300 MeV. To be identified as an e+, 
the track is required to exhibit at least two characteristics of positrons in heavy 
liquid. The signatures are: 

• bremsstrahlung with a 7 ^ e+e^ conversion; 

• spiralisation; 

• production of a (5 ray with energy comparable to the primary track; 

• annihilation with two 7 —>■ e^e~ conversions. 

After being scanned, measured and reconstructed, each event containing an 
e+ track is carefully analyzed, a fiducial volume cut applied and the search for 
at least one leaving negative {L~) track performed. The highest momentum L~ 
track in each event is interpreted as a muon leaving the chamber. Under this 
assumption, a kinematical analysis is applied in order to select genuine I'^CC 
interactions. 

2.3 The average di-lepton charm-production cross-section 

In order to combine all di-lepton data (see Tables ^ || and |^) it is necessary to 
make a bin by bin weighted average. Since different experiments have different 
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15 




40-60 
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30-60 
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60-100 


80 


0.71 ± 0.19 
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0.78 ± 0.24 


100-120 
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225 


0.95 ± 0.36 
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0.69 ± 0.13 








140-160 
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0.69 ± 0.13 
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E, 
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< E,, > 
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Oi-i+fo-cc 

(%) 


0-15 


7.5 


0.19 ± 0.07 


0-20 


10 


0.23 ± 0.11 


15-30 


22.5 


0.53 ± 0.09 


20-40 


30 


0.36 ± 0.11 


30-60 


45 


"•OO-0.08 


40-70 


55 


0.47 ± 0.14 


60-100 


80 


0.82ro.24 


70-100 


85 


0.61 ± 0.15 


1 on 1 f^n 


1 90 




1 nn 9nn 

lUU-ZUU 


ioO 


U.Do in U. 1 / 








200-300 


235 


0.53 ± 0.39 


CCFR 






Foudas et al. 
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E, 
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<E^> 
(GeV) 


oi-i+fo-cc 
(%) 




40 


0.54 ±0.05 


10-70 


50 


0.63 ±0.11 




70 


0.64 ± 0.04 


70-100 


80 


0.60 ± 0.07 




90 


0.66 ± 0.04 


100-150 


120 


'^-0.04 




110 


0.65 ± 0.03 


150-200 


170 


0.77 ±0.04 




140 


0.77 ± 0.04 


200-250 


220 


0.84 ±0.07 




180 


0.84 ± 0.03 


250-300 


270 






220 


0.86 ± 0.03 


300-350 


320 






270 


0.92 ±0.03 


350-400 


370 


o.97t;i:i^ 




330 


0.89 ± 0.05 


400-600 


450 


1.10 ±0.18 




450 


0.95 ±0.05 









Table 2: Di-lepton neutrino cross-section, normalized to the CC cross-section, 
as a function of the neutrino energy, for some (d to §, 0) electronic detector 
experiments. 
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<E^> 




(GeV) 


(GeV) 


(%) 


(GeV) 


(GeV) 


(%) 


0-25 


12.5 


0.20 ±0.04 


30-100 


70 


0.48 ±0.06 


25-50 


37.5 


0.51 ±0.08 


100-180 


150 


0.81 ±0.10 


50-100 


75 


0.60 ±0.10 


180-230 


200 


0.90 ±0.10 


100-200 


150 


0.60 ±0.15 









Table 3: Di-lepton neutrino cross-section, normalized to the CC cross-section, 
as a function of the neutrino energy, for some (0 to Q and ijlll) electronic 
detector experiments. 
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60-80 


70 
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140 
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30-100 


60 


0.38 ±0.13 


12-65 


50 


0.63 ±0.15 


100-230 


150 


0.85 ±0.27 


65-95 


80 


0.84 ±0.15 








95-145 


120 


0.72 ±0.12 








145-195 


170 


0.90 ±0.12 








195-290 


230 


0.84 ±0.12 








290-600 


350 


1.25 ±0.30 



Table 4: Di-lepton anti-neutrino cross-section, normalized to the CC cross- 
section, as a function of the neutrino energy, for the electronic detector experi- 
ments. 
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binning, we have to define weighting criteria to combine them all into an unique 
re-binning of the neutrino energy range: 0-600 GeV divided into 60 bins. The 
data-point in the j — th bin {{Xj,Yj),j e {1, . . . , 60}) is calculated once we have 
assigned a given weight to the data-point (xn- , ytk) of the A; — th bin of the i — th 
experiments, here-after called source bin, according to the following criteria: 

• the wider the source-bin {Axik ), the smaller the weight; 

• the wider the error {Ayik ), the smaller the weight; 

• the larger the fraction fijk of the j — th bin covered by the source bin, the 
larger the weight. 

It is then possible to make a bin by bin weighted average of the data-points 
by using the following formulae: 

• data-point weight 



Since CCFR data do not report the energy bins, it is not possible to combine 
them with the others at this stage of the analysis. 

Once the average has been computed, the CCFR data are combined with 
the other ones in the following way: each CCFR data-point is assigned to the 
bin containing its x-value, then a standard weighted average between the CCFR 
data-point and the previous combined data is calculated. In this case, the inverse 
of the Ay squared is used as weight. 

The average di-lepton single-charm production cross-sections for both ly and 
u are shown in Tables || and ||, and Figs. |l| and |2[ respectively. 

3 Charm-production studies with nuclear emul- 
sions 

3.1 Experimental issues 

So far only two experiments, E531 [|l2j and CHORUS have searched for 
inclusive charm-production through the direct identification of charm decays 
in the emulsions. The main advantage of these experiments is that, being the 




• data-point average 




'l2i,k ^ijk 



• data-point error 
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65 


0.53 ± 0.05 


205 


0.84 ± 0.05 






75 


0.61 ± 0.03 


215 


0.84 ± 0.05 
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0.60 ± 0.05 


225 


0.85 ± 0.03 






95 


0.64 ±0.03 


235 


0.81 ±0.06 






105 


0.66 ±0.07 


245 


0.84 ± 0.07 






115 


0.65 ± 0.03 


255 


0.86 ± 0.08 






125 


0.72 ±0.08 


265 


0.86 ±0.08 






135 


0.72 ± 0.08 


275 


0.91 ± 0.03 







Table 5: Average neutrino di-lepton cross-section, normalized to CC cross- 
section, as a function of the neutrino energy. 
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0.85 ±0.16 
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0.84 ±0.12 
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0.63 ±0.15 
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0.87 ±0.14 
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1.25 ±0.30 
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0.63 ±0.15 
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1.25 ±0.30 


55 


0.53 ±0.10 


195 


0.88 ±0.09 






65 


0.71 ±0.06 


205 


0.84 ±0.11 






75 


0.72 ±0.11 


215 


0.84 ±0.11 






85 


0.74 ±0.05 


225 


0.96 ± 0.08 






95 


0.74 ±0.13 


235 


0.84 ±0.11 






105 


0.90 ±0.07 


245 


0.84 ±0.12 






115 


0.77 ±0.13 


255 


0.84 ±0.12 






125 


0.80 ±0.14 


265 


0.91 ±0.08 






135 


0.86 ± 0.06 


275 


0.84 ±0.12 







Table 6: Average anti-neutrino di-lepton cross-section, normalized to CC cross- 
section, as a function of the neutrino energy. 
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Figure 1: Average neutrino di-lepton cross-section, normalized to CC cross- 
section, as a function of the neutrino energy 
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Figure 2: Average anti-neutrino di-lepton cross-section, normalized to CC 
section, as a function of the neutrino energy. 
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charmed particle identified through its decay, very loose kinematical cuts are ap- 
plied. This translates into a very good sensitivity to the slow-rescaling threshold 
behavior and consequently to the charm-quark mass. 

These experiments have a hybrid design. Electronic detectors locate the re- 
gions in emulsions where the neutrino interactions occur and contribute to the 
reconstruction of the event kinematics. Emulsions are used as active targets. 
They have the appropriate position resolution (less than 1 /xm) and granularity 
to detect short-lived particles through the visual observation of their decays. 
The main background for detection in emulsion comes from i^T" and A de- 
cays, and neutron, K'^ and A interactions without any visible nuclear break-up 
at the interaction point. For charged charmed- hadrons the main backgrounds 
are tt and K decays in flight, and the white kink (hadron interaction without 
any visible nuclear break-up) on any charged non-charmed hadron. Both the 
processes can be easily ruled out by applying a cut on the transverse momentum 
at the decay vertex {pr > 250 MeV). The contribution to the background of 
these processes is of the order of 

3.2 Emulsion data sample and selection criteria 

Although in the last years the automatic scanning of the nuclear emulsions had 
had an impressive development, the event search principle used in the former 
hybrid experiment E531 is still valid for present experiments. Once the neu- 
trino interaction is successfully reconstructed, tracks are searched for in the 
downstream end of the emulsion target by using the slopes and coordinates as 
predicted by the electronic detectors. If the predicted track is found in the most 
downstream emulsion sheet, it is then followed through the emulsion module up 
to the neutrino interaction vertex. 

Once a neutrino interaction is found, an attempt is made to find possible 
charmed-particle decays. In the former E531 experiment semi-automatic sys- 
tems were used, thus limiting the number of interactions to be studied. Nowa- 
days, very powerful automatic systems are available allowing for a larger statis- 
tics to be searched for. 

The final E531 data sample consists of 3855 neutrino interactions among 
which 122 charmed-particle decays are found The estimated background is 
0.2 and 3.6 events for neutral and charged charmed-particle decays, respectively. 

Using the weighted number of events, the relative charmed-particle produc- 
tion rate is: a{Vf,N c^i~ X)/(t{v^N n'X) = 4:.9tol%- Table [t] shows the 
energy dependence of this cross-section. The analysis of the CHORUS data is 
still under way. 

4 The world average charm-production cross- 
section 
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<E^> 


O'c/o'CC 


(GeV) 


(GeV) 


(%) 


0-10 


5 


'-'•'-0.6 


10-30 


20 


3.6 ±0.9 


30-60 


45 


4.9 ± 1.1 


60-150 


105 


10.1 ± 2.5 



Table 7: Inclusive neutrino charm-production cross-section, normalized to the 
CC cross-section, as a function of the neutrino energy. 



4.1 Prom the di-lepton to the total charm-production cross- 
section 

As discussed in Section |[ calorimetric and bubble chamber experiments cannot 
detect directly the charmed hadrons; instead only the lepton from semi-leptonic 
charm decay is measured. Without going into the details, we remind that the 
charmed hadron (C) cross-section can be connected to the charmed quark (c) 
cross-section via fragmentation functions, assuming factorization [p^ : 



da{i^N fj,-CX) _ d<7{vN ^ ^i-cX) 
dxdydzdp^ dxdy 

da{PN n+CX) _ d(j{vN fi+cX) 
dxdydzdp^ dxdy 



X 

h 



J2fkXD'^iz,p^^) 

h 

J2hxD^,i^,PT) (2) 



Here, D'^{z,p'^) is the probability distribution for the charmed-quark to frag- 
ment into a charmed hadron of type h{= D+, D+, A+) carrying a fraction z 
of the quark longitudinal momentum and transverse momentum pT with respect 
to the quark direction. The number fh is the mean multiplicity of the hadron h 
in neutrino charm-production. Analogously for anti-neutrinos, indicated by the 
barred quantities. Since only one c-quarks is produced in a CC interaction, one 
can set the normalization conditions as 

dz D'^{z,p^)dp'^ = dz 5|(z,p|)dp| = l (3) 
Jo Jo Jo Jo 



and 



The total cross-section for di-lepton charm-production, at a given energy, 
can then be written as 



a{iyN n-l+X){E) = a{vN ii-CX){E) x ^ fh{E) x BR^ 

h 
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Energy (GeV) 


/do 


/d+ 


fnt 


/a+ 


5-20 


0.32 ±0.11 


0.05 ± 0.06 


0.18 ±0.10 


0.44 ±0.12 


20 - 40 


0.50 ±0.08 


0.10 ±0.08 


0.22 ±0.08 


0.18 ±0.07 


40-80 


0.64 ± 0.08 


0.22 ±0.09 


0.09 ±0.08 


0.05 ± 0.04 


> 80 


0.60 ±0.11 


0.30 ±0.11 


0.00 ±0.06 


0.09 ±0.08 



Table 8: E531 prediction fraction results as obtained in Ref. 0. 



= aiiyN ^ l^rCX){E) X BRi{E) (4) 
where is the semi-leptonic branching fraction for the charmed hadron 

h. 

From the previous formulae it is clear that in order to go from the di-lepton 
charm-production to the inclusive charm-production, BR^ and fh{E) should be 
known. The semi-leptonic branching fractions are measured with an accuracy 
of about 10% jl^, while the /^(i?)'s have been measured directly only in one 
experiment, E531. However, a bias was found in the way E531 extracted fh{E). 
A detailed discussion on the re-fitting of these data with the bias removed is 
reported in Ref. |]l^ . For our purposes the measured energy dependence of fh 
is relevant, see Table |[ The error on BRi{E) has been computed as : 

(AS^,(£;))2 =^/2(^;).(ABi?0'+I]5fi'-(A/.(^;))'±2-^Bi?,.Bi?,.coz;,,,(£;^ 

i i j>i 

where covj^i{E) is the covariance matrix of the fh values, defined as 

cov,,j{E) = AME) ■ Afj{E) ■ con^j{E) 

and corij{E) is the correlation among the fh [ p^ . 

The average inclusive neutrino charm-production cross-section, normalized 
to the CC cross-section, as a function of the neutrino energy is shown in Table |^ 
and Fig. |[ 

4.2 The world average total charm-production cross-section 

The inclusive charm-production cross-section as measured in the E531 experi- 
ment and as derived from di-lepton data are shown in Fig. |^. The agreement 
between the two determinations is quite good, although the statistical errors 
are large. A better determination of the charmed fractions in the CHORUS 
experiment would reduce significantly the errors on the cross-section as derived 
from di-lepton data. 

In order to parametrise the inclusive charm-production cross-section as a 
function of the energy, we have performed a fit with a polynomial function 
(/(x')) requiring that the constant term is zero: 
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285 


9.05 ± 2.00 
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25 


r ro + i.4b 


165 


7.79 ± 1.64 


350 


9.68 ± 1.98 


35 


3.82+i:f^ 


175 


7.79 ± 1.64 


450 


10.11 ± 2.09 


45 


5.67 ± 1.06 


185 


8.63 ± 1.75 


550 


11.58 ±2.99 


55 


5.44 ± 1.12 


195 


8.11 ± 1.70 






65 


5.89 ± 1.19 


205 


8.84 ± 1.85 






75 


6.78 ± 1.25 


215 


8.84 ± 1.85 






85 


6.32 ±1.37 


225 


8.95 ± 1.82 






95 


6.74 ± 1.38 


235 


8.53 ±1.82 






105 


6.95 ±1.57 


245 


8.84 ± 1.92 






115 


6.84 ± 1.40 


255 


9.05 ± 2.00 






125 


7.58 ± 1.73 


265 


9.05 ± 2.00 






135 


7.58 ± 1.73 


275 


9.58 ±1.94 







Table 9: Average inclusive neutrino charm- production cross-section, normalized 
to the CC cross-section, as a function of the neutrino energy. 



ai ± Aoi 




ai 


a2 


aa 


a4 


as 


as 


0.23 ±0.03 


ai 


1.000 


-0.828 


0.266 


0.228 


0.028 


-0.204 


-(0.31 ±0.03) X 10"^ 


a2 


-0.828 


1.000 


-0.664 


-0.132 


0.096 


0.196 


(0.21 ±0.07) X 10~* 


as 


0.266 


-0.664 


1.000 


-0.463 


-0.168 


0.151 


-(0.71 ±0.02) X 10"^ 


a4 


0.228 


-0.132 


-0.463 


1.000 


-0.371 


-0.225 


(0.115 ±0.003) X 10"" 


as 


0.028 


0.096 


-0.168 


-0.371 


1.000 


-0.735 


-(0.71 ±0.05) X 10"" 


as 


-0.204 


0.196 


0.151 


-0.225 


-0.735 


1.000 



Table 10: Fitted parameters ai together with the correlation coefficients 
cor{ai,aj). 



f{x) = aix ± ± a^iX ± . . . + a„a;" 

Starting from rt = 6, the of the fit is very good, P{x^) > 0.99, and 
mildly dependent on n. Therefore, we have chosen the smallest power for n. 
We have studied the effect of the different parametrisations on the prediction of 
the inclusive charm-production cross-section and we have found results different 
of at most 5% from the ones obtained with n = 6. Therefore, we have assigned 
to each prediction a systematical error of 5%. The fitted parameters together 
with the correlation coefficients cor(ai,aj) are given in Table 10, The fitted 



function is also shown in Fig. |^. Notice that the present parametrisation holds 
in the range [0 300] GeV. 

As an application of our fit to the available data, we compute the expected 
D° production, normalised to the neutrino CC cross-section, in the CHORUS 
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experiment, whose neutrino beam has an average energy of 27 GeV. By con- 
voluting the neutrino spectrum with the cross-section shown in Fig. the 
expected average inclusive charm-production, normahzed to the neutrino CC 
cross-section, is (4.81 ± 0.16 \stat ±0.24 \syst)%, while the expected a{D^)/acc 
is (2.58 ± 0.26 \stat ±0.13 \syst)%- This has to be compared with the mea- 
sured value of 2.34 ± 0.15 ± 0.17% Notice that the former experiment 
E531, with an average neutrino energy of about 22 GeV, measured the ratio 

a{DO)/acc^2.19t'oiim- 

5 Conclusion 

Over the past 30 years, several experiments have accumulated data on the (anti- 
)neutrino induced single-charm production with different techniques: opposite- 
sign dimuons have been used by calorimetric experiments, di-leptons (/j,, e) by 
bubble chambers while only nuclear emulsion experiments have identified the 
charm through the visual observation of its decay. By combining all these data, 
a wide energy range (0-600 GeV) has been spanned, providing a good sensitivity 
to the threshold behavior of the cross-section. 

We have presented a statistical way of combining all the available data and hence 
the extraction of the world averaged single-charm production cross-section as 
a function of the neutrino energy. A comparison of the production cross- 
section, normalised to the neutrino CC cross-section, with the result recently 
presented by the CHORUS Collaboration shows a good agreement. 
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Figure 3: Average inclusive neutrino charm-production cross-section, normal- 
ized to the CC cross-section, as a function of the neutrino energy. The contin- 
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